EPJ manuscript No. 

(will be inserted by the editor) 



Quark Matter in QC 2 D 

Simon Hands 1 , Seyong Kim 2 and Jon-Ivar Skullerud 3 



o 
o 

(N 
Or 

in 



> 



Department of Physics, Swansea University, Singleton Park, Swansea SA2 8PP, U.K. 

2 Department of Physics, Sejong University, Gunja-Dong, Gwangjin-Gu, Seoul 143-747, South Korea 

3 School of Mathematics, Trinity College, Dublin 2, Ireland 

Received: date / Revised version: date 

Abstract. Results are presented from a numerical study of lattice QCD with gauge group SU(2) and two 
flavors of Wilson fermion at non-zero quark chemical potential (i > T. Studies of the equation of state, 
the superfluid condensate, and the Polyakov line all suggest that in addition to the low density phase of 
Bose-condensed diquark baryons, there is a deconfined phase at higher quark density in which quarks form 
a degenerate system, whose Fermi surface is only mildly disrupted by Cooper pair condensation. 
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1 Introduction 



The phase structure of QCD at large baryon density is 
one of the most fascinating areas of strong interaction 
physics, and yet a systematic calculational approach to 
this problem remains elusive. Lattice QCD simulation, the 
usual non-perturbative approach of choice, fails dismally 
because in Euclidean metric the quark action q~M(p)q, 
where M — If)[A] + /i7o + m with /i the quark chemical 
potential, results in a complex-valued path integral mea- 
sure detM when fi =/= 0. Since fi > promotes baryon 
current flow in the positive i-direction, the fundamental 
reason for this Sign Problem can be traced to the explicit 
breaking of time reversal symmetry. Because the measure 
no longer has an interpretation as a probability distribu- 
tion, Monte-Carlo importance sampling, the mainstay of 
lattice simulations, is completely ineffective in the ther- 
modynamic limit. 

It is instructive to ask what goes wrong when simula- 
tions are performed with a measure detM^M which is pos- 
itive definite by construction, as is the case for all practical 
fermion algorithms. It turns out that while M describes a 
color-triplet quark q £ 3, M* describes a conjugate quark 
q c 6 3. The model's spectrum thus contains gauge-singlet 
qq c states, indistinguishable from mesons at /x = 0, but 
carrying non-zero baryon number. As (i rises, baryonic 
matter first appears in the ground state (i.e. n q > 0) at an 
onset /i Q ~ ^TOtt, i.e. with an energy per quark comparable 
with the lightest baryon in the spectrum, which is degen- 
erate with the pion, rather than the physically expected 
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Only calculations performed with the 



correct measure det f M have cancellations among config- 
urations, due to the fluctuating phase of the determinant, 
which ensure that n q vanishes for \rn-n < fi < \m nuc i eon . 



For Two Color QCD (QC 2 D), ie. for gauge group SU(2), 
this bug is actually a feature. Since q and q live in equiva- 
lent representations of the color group, hadron multiplets 
contain both qq mesons and qq baryons. It is correspond- 
ingly straightforward to show that the quark determinant 
is positive definite for even Nf pQ. QC2D is thus the sim- 
plest model of dense strongly-interacting matter amenable 
to study with orthodox lattice techniques. Additionally, 
if there is a separation of scales <C m p in the spec- 
trum, then at low densities attention may be focussed on 
the Goldstone bosons of the system (both mesons and 
baryons) using chiral perturbation theory (%PT) The 
key result is that for fi > fi = a non- vanishing 

quark density n q > develops, along with a superfluid 
diquark condensate (qq) ^ 0. Just above onset, the sys- 
tem is thus a textbook Bose Einstein Condensate (BEC) 
formed from tightly bound scalar diquarks. 

Using the %PT prediction for n q (n) it is simple 
to develop the full equation of state, i.e. pressure p and 
energy density e q , at T = [3]: 



/4 



n q = 8N f f^(l-^j; 
p = J^n q d^ = AN f f^ + ^-2^) 
e q = -p + ^n q = 4N > ft U - ^ + 



(1) 



Here, / ff is a parameter of the model. Contrast this with 
another paradigm for cold dense matter, namely a degen- 
erate system of weakly-interacting massless quarks popu- 
lating a Fermi sphere up to some maximum momentum 
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(2) 



Superfluidity in this scenario arises from the condensation 
of quark Cooper pairs within a layer of thickness A centred 
on the Fermi surface, so that (qq) oc A\j? . 

Fig. n plots n q , p and e q from Q, each divided by 
the free field results (J5J, as functions of fi. On equat- 
ing pressures, this naive model, which ignores all non- 
Goldstone and gluonic degrees of freedom, predicts a first 
order deconfining transition from BEC to "quark matter" 
at fid f» 2.3/x with the choice f n — N c /6ir 2 . 



2 Simulation 

To test whether this prediction holds in a more systematic 
calculation we have performed simulations of SU(2) lattice 
gauge theory with Nf — 2 Wilson fermions with fi ^ [3] . 
The Wilson formulation is not obviously a stupid choice: 
Wilson fermions retain a conserved baryon charge; any 
problems with chiral symmetry should dominate in the 
\ow-k region of the quark dispersion curve, which lies at 
the bottom of the Fermi Sea and is hence inert; moreover, 
studies with free fermions show that saturation artifacts 
due to the complete filling of the first Brillouin zone ac- 
tually set in at higher values of fj, than is the case for 
staggered @]. Most importantly, the eigenvalue spectrum 
of the Wilson Dirac operator has the same symmetries as 
that of continuum QC2D. As shown in this fact per- 
mits an exact ergodic hybrid Monte Carlo algorithm for 
Nj = 2, with no requirement to take a fourth root, which 
may be problematic for /1 ^ [5]- The only novelty of our 
simulation is the inclusion of a diquark source term 

jqq = jK(-$i(x)C'y S T2iptr(x) + ipzix^^iix)) (3) 

in the dynamics, where subscripts label flavor and the 
Pauli matrix acts on color. As well as making the algo- 
rithm ergodic, setting j ^ mitigates the effect of IR fluc- 
tuations due to Goldstone modes in any superfluid phase, 



and of course enables direct estimation of the (qq) con- 
densate. 

Our initial study has been performed on an 8 3 x f6 
lattice using a standard Wilson gauge action, with param- 
eters (3 = 1.7, k — 0.178, and j = 0.04 (with a few points 
taken at j = 0.02, 0.06). Studies of the static quark poten- 
tial and the hadron spectrum at p, = yield a = 0.220fm, 
m^a = 0.79(1), and m^/mp = 0.80(1). 1 We thus expect 
the onset of baryonic matter at fj, a « 0.4. Thermody- 
namic observables are calculated as follows: quark density 
is given by a local operator via 



ainZ 
<9/i 



(4) 



As a component of a conserved current, it is immune from 
quantum corrections, but may be affected by artifacts due 
to a > 0, V < 00. We therefore prefer to quote our results 
in terms of n q /n l gg, where Tigg((i) is evaluated for free 
massless quarks on the same lattice. The pressure follows 
from an integral formula 
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(5) 



Note that although p is calculated purely in terms of quark 
observables, it is in principle the pressure of the system 
as a whole, although both continuum and thermodynamic 
limits must eventually be taken. Finally, quark energy den- 
sity is also estimated by a local operator 



£ q = ^: E (7o-lKlWVVf6-' , M7o + l)e ^^Vo^-O/ 

(6) 

this requires both subtraction of the \i = vacuum con- 
tribution, and a /i-independent but as yet unknown mul- 
tiplicative renormalisation. In what follows, therefore, the 
shape of the curve is in principle correct, but the overall 
scale still undetermined. 

Fig. [2] summarises our results. Both n q and p start 
to rise from zero at /za f» 0.3, although a careful j — > 



1 this corrects the value erroneously given in 
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extrapolation will be needed to pinpoint the onset with 
any precision. By [ia ~ 0.5 both quantities scale with 
\i in general accordance with free-field predictions, but 
with approximately twice the expected value. One expla- 
nation of this mismatch is that the system has formed a 

1 

Fermi sphere with /1 = Ep < kp oc n| , which could be 
attributed to a negative binding contribution to E from 
interactions. The quark energy density by contrast, in- 
creases more slowly than free-field expectations up to [ia w 
0.65, whereupon free-field scaling sets in rather abruptly. 
Another intriguing result [3] is that for 0.4 < [ia < 1.0 
the gluon energy density s g (identically zero in free-field 
theory) scales to quite high precision as /i 4 , the only phys- 
ically sensible possibility once fi/T 3> 1. Note that e g > 
entirely as a result of interactions with the background 
quark density, since this is the only means by which \i- 
dependence can arise. 

To elucidate what's happening, Fig. 3 plots both the 
superfluid order parameter (qq) divided by /i 2 , and the 
Polyakov line L. For /ia > 0.5 it is clear the system is 
in a superfluid phase, but what is remarkable is that at 
jua ks 0.6 there is a sudden transition to a regime where 
(qq) oc /1 2 , as expected for BCS pairing at a Fermi surface. 
At roughly the same point L rises from zero; although for 
theories with fundamental matter L is not strictly an or- 
der parameter, this is suggestive that at \ia « 0.65 there 
is a deconfining transition, beyond which the effective de- 
grees of freedom are best thought of as quarks (or even 
quasiquarks) , and not the scalar diquarks of xPT. 

3 Discussion 

Our initial study of thermodynamic quantities, and of 
the properties of the ground state, strongly suggests that 
QC2D at low temperature has at least two transitions as 
chemical potential fj, is raised. The first is between the 
vacuum and a phase of Bose-condensed tightly-bound di- 
quarks; the second, a relativistic analogue of the BEC/BCS 
crossover currently discussed in both strongly-correlated 
electron and cold atom systems, is a deconfining transition 
to a system of degenerate quarks, the Fermi surface being 
mildly disrupted by a Cooper pair condensate. Although 



QC2D clearly models nuclear matter unrealistically, its 
description of quark matter may well prove to have much 
in common with that of QCD. We are currently extend- 
ing our study to the hadron spectrum, and to finer lattice 
spacings to check that this conclusion is not due to lattice 
artifacts. Interesting results obtained from a study of the 
gluon propagator on the current system will be discussed 
elsewhere jHHH] 

Meanwhile it is hard to resist the temptation to spec- 
ulate on what a Two Color Star might look like. Fig. 
plots the energy per quark e q /n q versus \i using the data 
of Fig. 13 The most striking feature of this plot is the 
pronounced minimum at fia w 0.8, which is both robust 
(since it occurs even if corrections for a > 0, V < 00 are 
left out), and unexpected (since it does not occur for the 
model EoS of Fig. We infer that any large object as- 
sembled from a fixed number of QC2D quarks, such as a 
star, will have the bulk of its interior in the neighbour- 
hood of this minimum, which as Fig. 3 shows, means that 
the object would in effect be a quark star formed from de- 
confined matter. Somewhat speculatively, we have labelled 
the different regions of the //-axis with the corresponding 
layers of the star, although a quantitative solution for the 
radial profile must await correctly-normalised calculations 
of the energy densities e q and e g . 



References 

1. S. Hands, I. Montvay, S. Morrison, M. Oevers, L. Scorzato 
and J.I. Skullerud, Eur. Phys. J. C 17 (2000) 285. 

2. J.B. Kogut, M.A. Stephanov, D. Toublan, J.J.M. Ver- 
baarschot and A. Zhitnitsky, Nucl. Phys. B 582 (2000) 477. 

3. S. Hands, S. Kim and J.I. Skullerud, prXiv:he p-lat /0604004 
Eur. Phys. J. C (in press). 

4. W. Bietenholz and U.J. Wiese, Phys. Lett. B 426 (1998) 
114. 

5 . M. Golterman, Y. Shamir and B. Svetitsky, 
|arXiv:hep -lat /0602026 

6. J.I. Skullerud, talk at this conference. 



